Multi-Domain Spectral Method for Initial Data of 
Arbitrary Binaries in General Relativity 



Marcus Ansorg 

Max-Planck-Institut fur Gravitationsphysik, Albert-Einstein-Institut, Am 
Muhlenberg 1, D-14476 Golm, Germany 

E-mail: marcus.ansorg@aei.mpg.de 

Abstract. We present a multi-domain spectral method to compute initial data of 
binary systems in General Relativity. By utilizing adapted conformal coordinates, the 
vacuum region exterior to the gravitational sources is divided up into two subdomains 
within which the spectral expansion of the field quantities is carried out. If a component 
of the binary is a neutron star, a further subdomain covering the star's interior is added. 
As such, the method can be used to construct arbitrary initial data corresponding 
to binary black holes, binary neutron stars or mixed binaries. In particular, it is 
possible to describe a black hole component by the puncture ansatz as well as through 
an excision technique. First examples are given for binary black hole excision data 
that fulfill the requirements of the quasi-stationary framework, which combines the 
Conformal Thin Sandwich formulation of the constraint equations with the Isolated 
Horizon conditions for black holes in quasi-equilibrium. These numerical solutions were 
obtained to extremely high accuracy with moderate computational effort. Moreover, 
the method proves to be applicable even when tending toward limiting cases such as 
large mass ratios of the binary components. 



PACS numbers: 04.25.Dm, 04.20.Ex, 04.70.Bw 



1. Introduction 



The past year has seen considerable progress in the numerical calculation of the 
dynamics of binary black holes. Different groups were able to demonstrate computations 
describing multiple orbits, the merger as well as the final ring-down phase with the 
corresponding gravitational radiation emitted PElEllligillTJIElEani^ 
While most of these groups use initial data arising from the puncture ansatz [15J, 
different initial data, which result from the quasi- stationary framework of Conformal 
Thin Sandwich equations [HI |T7] (see also review [18]) and Isolated Horizon boundary 
conditions [19, 20J (see also [21]; for reviews on the Isolated Horizon formalism see 
[22 | [23 | 121] ) have also been chosen. As a consequence, the numerics of dynamical binary 
black holes has reached a stage in which important astrophysical implications can be 
addressed, e.g. the accurate calculation of the final black hole's "kick" velocity which 
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appears for unequal-mass binaries due to the asymmetric loss of linear momentum in 
the gravitational radiation emitted during the merger (see e.g. [T3]). 

The initial data used in these calculations do not represent an exact astrophysical 
initial situation of distant orbiting black holes, but only an approximation to this. 
In particular, through a specific (often mathematically convenient) fixing of "free" 
components in the initial data construction, an undesired and uncertain amount of 
artificial gravitational radiation is contained in the data. In the calculation of the 
dynamics of neutron stars as well as that of mixed binaries (consisting of a black hole and 
a neutron star) the same issues arise. It is therefore necessary to check the numerically 
calculated astrophysical results with respect to a large variety of different types of binary 
initial data. 

Indeed, a large variety is available because of the fact that, in order to construct a 
specific initial data set, one has to 

(i) choose a formulation of the constraint equations, 

(ii) fix the free components in the data, 

(iii) for black holes: choose between the puncturdjl and the excision method^], 

(iv) for black holes with excision technique: fix the formulation of the boundary 
conditions needed at the excision surface. 

In this paper we introduce a new multi-domain spectral method to compute 
arbitrary initial data of binary systems in General Relativity. Apart from the freedom 
to provide any type of data, we aim at an extremely high numerical accuracy close to 
machine precision, even when tending toward limiting cases such as large mass ratios or 
large distances of the binary components. Because of their exponential convergence 
properties, we utilize pseudo-spectral methods, which also have been used e.g. in 
[251 12H1 [271 12HI [291 [2H [20l EQl EH [321 [33l EH [35] for the construction of initial data. 

Conformal coordinate mappings for the construction of initial data in General 
Relativity have been explored in the literature, see e.g. [361 EH EH El]- The basic 
feature of the method described here is a specific conformal coordinate mapping which 
permits the identification of the entire spatial domain as the image of a small number of 
compact subdomains, introduced such that the data are analytic there. Consequently, 
a rapidly converging spectral expansion ariseijj. In particular, the exterior vacuum 
region is divided up into two subdomains, and for each neutron star component of 
the binary a subdomain covering all of the star's interior is added. Thus we obtain a 

\ In the puncture methods, a special pole- like structure of the singularity (called the "puncture" ) inside 
the black hole is assumed which can be taken into account by a specific ansatz for the initial data. 
Therefore the relevant space for the constraint equations is all of M 3 . 

§ The excision techniques solve the constraints only in the exterior of an excised spheroid within which 
the singularity is located. 

|| Note that the corresponding coordinate curves are in general not smooth at common domain 
boundaries, see figures [3] -[8] This does not lead to a loss of accuracy in the numerical scheme since the 
spectral expansions are performed separately in each subdomain and the communication between them 
is realized through transition conditions (continuous behaviour of functions and normal derivatives). 
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two- domain spectral method for binary black holes, a three-domain spectral method 
for mixed binaries and a four-domain spectral method for neutron star binaries. The 
method can be used to compute the following types of initial data: 

(i) Binary black hole puncture data 

(ii) Mixed binary black hole puncture-excision data 

(iii) Binary black hole excision data 

(iv) neutron star - black hole binary data; black hole: puncture 

(v) neutron star - black hole binary data; black hole: excision 

(vi) Binary neutron star data 

The paper is organized as follows. In section [2] the conformal mapping in question 
is introduced. Section [3] describes in which manner the partition into subdomains as 
well as their individual mappings can be undertaken for each one of the above types 
of initial data. The pseudo-spectral scheme which takes all subdomains into account is 
explained in section HI Finally, examples are given in section [5] for binary black hole 
excision data that fulfill the requirements of the above quasi-stationary framework. The 
corresponding code is of low computational cost, allowing one to do the calculations on 
a laptop computer. 

Note that we use units in which the speed of light as well as the gravitational 
constant are equal to one. An overbar denotes complex conjugation. 

2. The conformal coordinate mapping 

Our starting point is a Cartesian-like coordinate system covering all of the initial 
t = constant hypersurface. We denote these coordinates by (x, y, z) defined on the 
interval (—00,00). The gravitational sources are assumed to be located at finite 
coordinate values, thus describing an asymptotically flat spacetime. In particular, we 
choose the coordinates such that the binary source be aligned along the x-axis and the 
y-z-plane be located between the binary components. That is, for binary black hole 
puncture data, we place the punctures at the coordinate values 

x = ±b, y = = Z , (6>0). (1) 

In the case of sources that are either black holes, which are given by horizon boundary 
values at the surface of an excised spheroidal shell, or spheroidal neutron stars, we 
describe the surface of the spheroidal shells parametrically by 

x = S±($, f) cost? ± b, 

y = S±(d, p) sin $ cos (2) 

z = S±(d, p>) sin $ sin y), 

where the positive smooth functions S± can be expanded with respect to spherical 
harmonics Y^ifi, tp). For the examples considered in this paper we will refer to the 
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C = i,c = -b 




-b x b c = o, X C = l, 



c = oo c = b 

Figure 1. Illustration of the conformal mapping c = ife (C 2 + C~ 2 ) with c = x — \p 
and C = X + i R. Left and right panel show a (ip = constant)-slice in cylindrical 
coordinates (x, p) and the corresponding domain in the coordinates (X, R) respectively. 
The locations of the punctures are marked with crosses, and the curves around these 
points indicate spheroidal shells for excised black holes or neutron stars. 



simplest case, in which the shells are exactly spherical. However, the method also 
allows more sophisticated as well as unknown surface functions which, similar to the 
"free" boundary value problems for axisymmetric and stationary fluid configurations, 
are determined through some quasi-stationary condition imposed on the surface of the 
shellaS. 

If we are dealing with mixed binaries consisting of a puncture black hole and either 
a further excised black hole or a neutron star, we place the puncture at one of the points 
(CQ) and the shell according to (J2J) on the other side of the x-axis. The corresponding 
setting is illustrated in the left panel of figure [Q 

The introduction of appropriate conformal coordinate mappings has been explored 
in [30] and [31] in order to calculate binary black hole puncture and excision data 
respectively. For punctures the key idea was a "folding" of the x-axis about their 
locations ([I]), in order to obtain coordinates in which the data are analytic there. In 
a similar fashion, the "folding of infinity" in |3T] yields binary black hole excision data 
that are regular at the location of compactified infinity. In this paper we combine the 
two concepts and propose a conformal coordinate mapping that realizes the following 
requirements: 

(i) Compactification of spatial infinity, 

(ii) Regularity of the data at compactified infinity, 



% For neutron stars, the vanishing of the energy density at the surface provides such a condition, see 
[30] for a review. 
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(iii) Regularity of the data at the points (PQ) (for puncture and neutron star data), 

(iv) Identification of subdomains within which the data are analytic, 

(v) Mapping of these subdomains onto rectangular boxes within which the spectral 
coordinates 

(A,B,tp) e [0,1] x [0,1] x [0,2*) (3) 

are defined. 

We take cylindrical coordinates (x, p, tp) about the x-axis, i.e. 

x = x, y = pcos(p, z = psintp, (4) 
in order to form the complex combination 

c = x — ip. (5) 
The above requirements are met by the conformal transformation 



where 



Writing 



C = X + iR. 

explicitely gives 
1 



(6) 
(7) 



x 



{R 2 + x 2 y 
1 

{R 2 + x 2 y 
1 



X 2 -R 2 
2 ' 

RX cos(f, 
RX sin ip. 



(8) 



.(R 2 + x 2 Y 

The conformal mapping ([6]) is illustrated in figure [TJ It yields the entire spatial domain 
as the image of the cross-product of a quarter unit circle (on which X and R are defined) 
and the interval [0, 2n) over which the angle (p varies, i.e. 

c : {(X, R) : X > 0, R > 0, X 2 + R 2 < 1} x [0, 2vr 

In particular, 



U oo. 



(9) 



C 
C 

c 
c 




i 
1 



rV^(l + i) 



oo 
-b 
+b 




(10) 



The axes X = and R = are mapped to the sections x < — b and x > b of the x-axis 
respectively. The section between the punctures, — b < x < b (p = 0), is the image of 
the quarter circle curve. For X = R we obtain the p-axis. Figure Q] illustrates how the 
mapping folds the x-axis about the punctures as well as about infinity. 
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Figure 2. Illustration of the spectral mapping ([T3|) - (fT6|) for single domain binary 
black hole punctures. The circles X 2 + R 2 — constant correspond to constant A- 
coordinate lines (X 2 + R 2 — e^ = (1 — A) /{I + A)) and the rays X/R = constant are 
described by constant i?-values (X/R = cot(?y/2) = 1/B — 1). The right panel shows 
the corresponding curves in the cylindrical coordinates (x,p). 



For the regularity issues at compactified infinity consider the inverse radius 
1 1 1 2CC 



r \c\ b\l + C 4 |' 

which is regular at C = 0. Likewise, consider the distances 

r± = v^T^ + P 2 = \cTb\ = t~ |l =F C 2 \ 2 (12) 

in the vicinity of the puncture points (JTJ). They are regular at C 2 = ±1 which ensures 
the regularity of puncture data there. 

In cases with excised black holes or neutron stars, the spheroidal shells are indicated 
in figure [1] by the curves around the puncture points. Corresponding to the surface 
functions S±, these curves may depend on ip. The figure provides an intuitive partition 
of the spatial domain, dividing up the exterior vacuum region into two subdomains and 
taking for each neutron star component of the binary a subdomain covering all of the 
star's interior. The specific mappings of the spectral coordinates (A, B, tp) onto these 
subdomains for the various cases is the subject of section [31 

3. Spectral mappings for binary systems 

3.1. A single domain for binary black hole punctures 

For binary black hole puncture data, the quarter unit circle can be covered with a single 
coordinate patch using polar coordinates, see figure [2j In fact, with the conformal 
transformation 

C = e c/2 , (13) 
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Figure 3. Illustration of the spectral mapping (fl7|) - (fT9|) for double domain binary 
black hole punctures. In the vicinity of the punctures, the coordinates A and B 
resemble radial and angular coordinates respectively. The corresponding curves of 
constant A- and B-values are shown in the left and right panel. Specifically, in this 
example the coordinate lines A = constant form a dense mesh in the vicinity of the 
puncture x = —b, in order to resolve steep gradients. 



where 



we obtain 



C = £ + % (14) 



c = b cosh C, (15) 

i.e. we recover the mapping used in [SU] for the calculation of binary puncture data in 
a single spatial domain. Note that here £ < whereas £ > in [30], according to the 
fact that in that paper the exterior of a half unit circle has been mapped onto M. 3 . 

Hence we may write the coordinates £ and rj in terms of spectral coordinates in a 
similar fashion: 

£ = — 2artanh A , 

(16) 

r) = f + 2arctan(2 J B - 1) (A, Be [0,1]). 

The single domain spectral method for binary black hole puncture data has been used 
by several groups for the dynamical evolution of binary black hole space-times. For 
comparable masses and moderate distances of the black holes, this method yields a 
good convergence rate of the solution as one increases the numerical resolution. 



3.2. Two domains for binary black hole punctures 

If one considers large distances and/or large mass ratios of the black holes, the gradients 
of the gravitational field quantites grow large in the vicinity of the punctures. As 
a consequence, the spectral convergence rate of the single domain method drops 
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significantly. In order to retain the convergence rate in these coordinate mapping 

can be introduced where one of the coordinates resembles a measure of the coordinate 
distance from the puncture within the high gradient regime. The rate by which the 
distance scales with this coordinate can be adjusted such that an appropriate resolution 
of the steep gradients is achieved. In this manner a double domain method for binary 
black hole punctures arises. 

A possible transformation which maps the spectral coordinates (A, B) e [0, l] 2 onto 
the upper vacuum section of the quarter unit circle is given by (see figure [3]) 



C 



Bexp ( ) + (1 - B)(l - A_ 



4 

Likewise, the lower section is obtained by 



;i7) 



C = Be^(jA + )+(l-B)(l-A + ). (18) 

In these formulas the rescaled spectral coordinates A± appear which provide the desired 
scaling of the coordinate A with the coordinate distance from the puncture. They 
depend on b, the half coordinate distance of the punctures, and some quantity m± that 
represents a local mass (e.g. the puncture's bare mass [15] or an appropriate horizon 
mass HHSgj+l 

sinh(Aloge±) m± , , 

A± = . , r-, e± - — . (.lyj 

sinn (log e± J o 

Similar logarithmic rescalings to resolve steep gradients have been explored in |43j. 
3.3. Two domains for excised binary black holes 

In the case of black hole excision data, the transformations (|T7|) . (fT8j) need to be modified 
in order to take the spheroidal excision surfaces into account. As in [31] we utilize 
bispherical coordinates [HI H5] 

C = 4 artanh (C 2 ) , (20) 

in which 

3fc(C) = constant ^0, < 3(C) < vr 

corresponds to exact spherical surfaces about the points ([I]). Allowing general shapes, 
we may describe the surfaces by 

Cl{B,ip) = tanh 



r r ±o- ± (B,p)>0, (21) 

where the functions a± are related to S± introduced in (j2J) and take care of the excised 
spheroidal surface shape, which either can be prescribed or which may be unknown and 
determined through a quasi-stationary condition to be imposed. 

+ Here the subscripts ± refers to the puncture at the coordinate location x = ±b(p = 0). 
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Figure 4. Illustration of the spectral mapping for double domain mixed binary 
black hole puncture-excision data. The mapping is performed according to the 
transformations (fl7|) and l|22p in the upper and lower part of the quarter unit circle 
respectively. In this example, the excision surface is a sphere. Left and right panel 
display the coordinate curves of constant A- and f?-values. 




Figure 5. Illustration of the spectral mapping for double domain binary black hole 
excision data. The mapping is performed according to the transformations (|22[) in the 
upper and lower parts of the quarter unit circle. In this example, the excision surfaces 
are spheres. Left and right panel display the coordinate curves of constant A- and 
B- values. 
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Figure 6. Illustration of the three-domain spectral mapping for mixed binary black 
hole-neutron star data, where the black hole is described by the puncture method. 
The mapping is performed according to the transformations ([23]) , (|22[) and ([18)) in the 
upper and lower part of the quarter unit circle respectively In this example, the star 
surface is a sphere. Left and right panel display the coordinate curves of constant A- 
and B-values. 



In terms of these functions the mapping of the upper and lower vacuum sections of 
the quarter unit circle can be written in the following form: 

C= (l-A±)[C±(B,<p)-BC±(l,<p)] 



+ .Bexp 



[C±(B,<p) 
i (^A ± + (1 - At) 



argC±(l,v?) 



(22) 



Again the rescaled spectral coordinates A± appear to render steep potential gradients, 
see flUD. 

The puncture mappings (TTTj) , (|T8|) can be used together with the excision mappings 
( I2ip . ( 1221) . in order to calculate mixed binary black hole puncture-excision data, see 
figure HI Double excision black hole binaries are obtained by mapping both subdomains 
according to (12T1) and (122)) , see figure [51 



3.4- Domains for neutron stars 

For the calculation of neutron star data we may take (122)) in order to describe the 
surrounding exterior vacuum domain. The inner region adds another subdomain. In 
analogy to expression (122)) . we find a corresponding mapping, through which inner 
regions (with the same surface shape (T2T)) ) are written in terms of the spectral 
coordinates: 



C= (1-A_)[C-(B^)-BC_(1 )V )] 

+ 5exp [i(|A_ + (l-A_)argC_(l, </>))] + iA_(l-5) (23) 
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Figure 7. Illustration of the three-domain spectral mapping for mixed binary black 
hole-neutron star data, where the black hole is described by the excision method. The 
mapping is performed according to the transformations ([23]) and ([22]) in the upper and 
lower parts of the quarter unit circle. In this example, the star and excision surfaces 
are spheres. Left and right panel display the coordinate curves of constant A- and 
B- values. 

and 

C= (1 - A+) [C + (B,tp) - B C+(l,<p)] 

+ £?exp [i (1 - A+) arg C+(l, <p)] + A+(l- B) (24) 

for stars at the negative and positive sections of the x-axis respectively. In figures [6] and 
[7] the mapping for mixed binary data, consisting of a black hole and a neutron star, is 
illustrated; figure M displays it for binary neutron stars. 

4. The numerical scheme 

The numerical scheme for calculating binary initial data in the various cases generalizes 
the methods presented in [31] for binary black hole excision data. Most of the techniques 
are similar, and we refer to that paper for the details. Differences arise in the case of 
neutron star data where three or four domains are involved. Moreover, conditions with 
respect to the surface shapes as well as certain required physical parameter relations 
need to be taken into account. 

In our spectral approximation we write all potentials, that are involved in the elliptic 
boundary value problem, in terms of Chebyshev expansions with respect to A and B 
and Fourier expansions with respect to ip. For the spectral gridpoints we choose the 
extrema of the Chebyshev polynomials so as to have gridpoints lying on the boundaries: 
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Figure 8. Illustration of the four-domain spectral mapping for binary neutron star 
data. The mapping is performed according to the transformations (|23|) . ([22]) and 
([24]) . (|22|) in the upper and lower part of the quarter unit circle respectively. In this 
example, the star surfaces are spheres. Left and right panel display the coordinate 
curves of constant A- and i?-values. 



where 



Bk 



<Pi 



sin 
2ttZ 



2(n B - 1) 



(25) 



0<j<n A , 0<k<n B , < I < n v . (26) 

The numbers nA,n B and describe the spectral expansion orders of our scheme. Here 
ua may be different in the various subdomains, in order to allow higher resolution in 
regimes with steep gradients. However, n B and n v are equal in all subdomains, in order 
to have the same spatial grid points on both sides of the domain boundaries. 

The potential values at the above gridpoints are combined in a vector describing the 
spectral approximation of the desired solution. Moreover, we include in this vector the 
surface shape values a±(B k , (pi), if they are to be determined. Likewise, we add certain 
parameters of the solution that are constrained through specific parameter relations. As 
an example of this, we may take the orbital angular velocity Q of the binary system as 
a further unknown and determine it through the equality of Komar- and ADM-mass, a 
condition that arises within a quasi-stationary framework. 

Having set up this vector, the collection of 

(i) elliptic equations valid within the subdomains, 

(ii) transition conditions to be imposed at common borders of the subdomains, 

(iii) specific regularity conditions to be fulfilled at the x-axis, 

(iv) specific regularity conditions at C 2 = ±1 (for puncture and neutron star data), 
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Figure 9. The convergence of the ADM- mass corresponding to binary black hole 
initial data sets with vanishing orbital angular velocity fl — 0. The geometrical 
parameters of the individual configurations are given by b = 5g_ and g+/ Q- G 
{1; 10; 100; 1000} where g± describes the coordinate radius of the excised spherical 
shells. For these axisymmctric calculations tia = tib = n and n v = 3 have been 
chosen (in both subdomains). We compared the corresponding results for the ADM- 
masses to those of reference solutions with n 6 {50; 60; 70; 80} for the above choices of 
g+/g-. 
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(v) specific asymptotic boundary conditions at C = (i.e. at spatial infinity), 

(vi) excision surface boundary conditions (for black hole excision data), 

(vii) surface shape conditions (for unknown, quasi- stationary black hole or neutron star 
boundaries) , 

(viii) parameter relations (for determining specific physical parameters of the solution), 

yields a discrete non-linear systemQ which is solved by Newton-Raphson iterations. 
As in [3Ul EI] the linear step inside this solver is performed with the preconditioned 
'Biconjugate Gradient Stabilized (Bi-CGSTAB)' method [16]. A good convergence of 
this method requires a so-called preconditioning, which we construct in complete analogy 
to [3TH I3"T] through a second order finite difference representation of the Jacobian matrix 
of the non-linear system. In particular, the preconditioner consists of successive plane 
relaxations with respect to the planes tp = constant = <pi. 

For the convergence of the Newton-Raphson scheme a nearby initial 'guess' is 
essential. It can be taken from previously calculated numerical as well as from known 
analytical solutions. Typically, only a moderate number of iterations within both the 
Newton-Raphson and the Bi-CGSTAB methods are needed to calculate the desired finite 
spectral approximation of the solution. 
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5. First examples on binary black hole excision data 

The representative calculations of binary black hole excision data in [31] are specific 
examples of the general methods presented in this paper. These initial data fulfill the 
requirements of the quasi-stationary framework, which combines the Conformal Thin- 
Sandwich Decomposition with the Isolated Horizon framework. 

In the Conformal Thin-Sandwich Decomposition the evolution of the metric 
between two neighbouring slices t = constant is considered. For quasi-equilibrium 
binary data, the constraint and evolution equations of General relativity are treated 
in a corotating frame of reference in which the time derivatives of the metric quantities 
are assumed to be small initially. The corresponding Hamiltonian and momentum 
constraints are elliptic equations to determine the initial data, i.e. the 3-metric and the 
extrinsic curvature of the slice. An additional elliptic equation follows from prescribing 
some value for the time derivative of the trace of the extrinsic curvature. It can be used 
to provide a specific initial time-slicing of the data. 

This quasi-stationary formulation is completed by a set of boundary conditions 
that control the data at the excision boundaries and at infinity. The excision boundary 
conditions are given through the Isolated Horizon framework and describe black holes in 
a quasi-equilibrium state, i. e. as instantaneous non-expanding horizons. In particular 
the following is required: 

(i) Within the initial slice, each excision boundary is an apparent horizon, i.e. a two- 
dimensional hypersurface with S 2 topology and the property that the outgoing null 
vectors k possess vanishing expansion. 

(ii) Initially, the apparent horizon is tracked along k and its coordinate location does 
not move in the time evolution of the data. 

The examples presented in this paper correspond to so-called conformally flat and 
maximal Conformal Thin Sandwich data sets of binary corotational black holes that obey 
the requirements of the Isolated Horizon framework. These data have been previously 
calculated in [201 G3], see also [21] for further details. 

An ingredient that enters this formulation is the orbital angular velocity Q of 
the binary system. At first we consider Q = which results in specific generalized 
Misner-Lindquist initial data (they do not represent two orbiting black holes in a quasi- 
stationary state). As an indicator for the accuracy achieved we monitor the convergence 
rate of the ADM- mass Madm- We find it to be geometric, as exhibited by a roughly 
linear decrease of the error in figure [H Note that almost machine accuracy is reached for 
all ratios considered, thus proving that the method is well suited to the case of extreme 
mass ratios (and likewise to situations in which the coordinate distance parameter b is 
extremely large). 

The second example corresponds to initial data of two corotational black holes 
in a quasi- stationary orbit. The orbital angular velocity Q is obtained by requiring 

* The equations and conditions are considered at the discrete gridpoints (|25|) . 
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Figure 10. The convergence of the ADM-mass corresponding to a corotating binary 
black hole initial data set in a quasi-stationary orbit. The geometrical parameters 
of the configuration are given by b — 5g_ and g + = Q- where g± describes the 
coordinate radius of the excised spherical shells. The data are characterized by the 
equality Madm = Mk through which an orbital angular velocity of f2 ss 0.036975/g_ 
emerges. For these calculations the spectral resolutions ua = Ub = 2n v + 1 = n have 
been chosen (in both subdomains). We compared the corresponding results for the 
ADM-mass to that of the reference solution with n = 51. 




the equality of the ADM-mass Madm and the Komar mass Mk which can be defined 
by a suitable surface integral at infinity (see [27J, EEJ HSJ EUJ EI]). Again a geometric 
convergence rate of the resulting ADM-mass arises, see figure ITUl 
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